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ABSTRACT 


Non-Destructive Evaluation (NDE) involves the characterization of material properties 
and/or defects without detrimental effects on the structure examined. NDE can be 
performed using ultrasound, x-rays, microwaves or magnetic flux. Ultrasound has a wide 
range of applicability in performing NDE of matenals and structures. In Ultrasonic Non- 
Destructive Testing and Evaluation, use is made of the properties of elastic waves 
propagation in solids in order to detect the defects and materials inhomogeneties. In 
Ultrasonic NDE, the different features of ultrasonic waves, namely frequency, amplitude, 
time of flight etc can be used for detection and sizing of internal defects. Ultrasonic 
Tomography (UT) is used to characterize a defect or inhomogeneity by Time of Flight 
(TOF) information. The accuracy of reconstruction depends on the correctness of TOF 
determination and also on the specific algorithm used. 

The present work is on the use of TOF measurement and subsequent image 
reconstruction. In this study, the thermo-elastic waves induced in glass/epoxy composites 
with metal inserts by Laser Based Ultrasonics (LBU) are considered. The Time of Fli^t 
has been used for tomographic reconstruction of a specimen having defect using 
Algebraic Reconstruction Technique (ART) and Multiplicative Reconstruction 
Technique (MART) algorithms. The Tomographic reconstructions using different values 
of relaxation parameter and initial guess have been used to discuss the nature and 
suitability of these algorithms for defect reconstruction. Finally, Simulation study is 
conducted to check the performance of various algorithms and review the experimental 
results. 
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Chapter 1 

INTRODUCTION 


1.1 INTRODUCTION 

The present work on the identification and reconstruction of defects in composite 
matenals requires a thorough understanding of behavior of composite matenals, 
ultrasonic wave propagation and analysis from expenmental and theoretical viewpoint 
A bnef introduction to related topics is given below 

Demands imposed by modem advanced technologies on matenals have become so 
diverse and severe that single component conventional matenals often cannot meet them 
It IS necessary to combine two or more matenals to make a composite to meet the 
demands Constituents not only contnbute their respective share, but their action 
transcends the sum of the individual properties and provides new performance 
unattainable by constituents acting alone Composites can be formally defined as 
matenals having two or more chemically distinct constituent matenals or phases with 
clear interfaces on a macro scale between them 

The fiber reinforced matenals are probably the most important class of composite 
materials which are widely used m structural applications particularly when weight 
saving IS crucial These matenals offer high strength-to-weight and modulus-to-weight 
ratios apart firom their properties of excellent corrosion resistance, formability and 
mechanical properties, which can be controlled to a significant degree 

Non-destructive testing (NDT) may be defined as the science of examimng objects, 
matenals or systems, m order to determine their fitness for certain purposes, without 
impainng their future usefulness or their desirable properties The term is generally 
applied to non-medical investigations of matenal integnty Non-destructive evaluation 
(NDE) of matenal plays a vital role in quality assurance of structural members dunng 



manufactuimg stage and operating life In general, the vanous NDE techniques can be 
placed into two categones active and passive The active techniques are those where 
something is introduced into or onto the specimen and a response is expected if a defect 
IS present where as passive techniques are those that monitor or observe the item during 
either a typical load environment or a proof cycle and attempt to determine the presence 
of a defect through reaction of specimen Magnetics, ultrasonics, radiography fall m the 
category of active NDE and acoustic emission, noise analysis, leak testing etc are the 
examples of passive NDE Being one of the most commonly used non-destructive testing 
(NDT) methods, ultrasonic testing (UT) is rapidly developing m recent years 

Ultrasonic inspection involves impinging a low energy, high frequency stress pulse mto 
the matenal under inspection and examimng the subsequent propagation of this energy It 
not only detects bulk and surface flaws, but also obtains mformation on matenal 
microstructure to determine engineenng properties, such as elastic moduli and ultimate 
tensile strength However, traditional ultrasound requires liquid or contact coupling for its 
transmission, making it difficult or impossible to apply m many mdustnal situations This 
occurs m particular on curved parts and on parts at elevated temperature, a situation 
widely found in mdustnal products and dunng the processing of mdustnal matenals 

In conventional method, a piezoelectnc transducer (or probe) generates ultrasonic waves, 
which propagate in the elastic medium and are detected either by the same (pulse-echo) 
or by a different transducer (through transmission) Although piezoelectnc transducers 
are commonly used for non-destmctive testing, several problems are associated with the 
requirement that they have to be bonded to the test matenal with an acoustical impedance 
matching coupling medium For velocity measurements, which are necessary for matenal 
thickness measurements and to locate the depth of defects, the coupling medium can 
cause transit time errors Due to the partial transmission and partial reflection of the 
ultrasonic energy m the couplant layers, there may be a change of shape of the waveform, 
which can further affect the accuracy of the velocity measurement This can also lead to 
senous errors in absolute attenuation measurements, which is the reason that so few 
reliable absolute measurements of attenuation are reported m the scientific literature 
Therefore, a method of non-contact generation and detection of ultrasound is of great 



practical importance Several such techniques are presently available in vanous stages of 
development, namely capacitive pick-ups, electromagnetic acoustic transducers 
(EMATS), laser beam optical generators and detectors, and more recently air (or gas) - 
coupled ultrasonic systems However, as the name implies, capacitive pick-ups cannot be 
used as ultrasonic generators and, even when used as detectors, the air gap required 
between the pick-up and test structure surface is extremely small, which m essence 
causes the device to be very nearly a contact one The problem with EMATS is that the 
efficiency of ultrasound generation and detection rapidly decreases with lift-off distance 
between the EMATS face and the surface of the test object They can obviously be used 
only for examination of electncally conductive matenals Because of the physical 
processes involved they are much better detectors than generators of ultrasound Laser- 
based ultrasound (LBU) generation and detection overcomes all of these problems and 
make truly non-contact ultrasomc measurements in both electncally conductive and non- 
conductive matenals, m matenals at elevated temperatures, in corrosive and other hostile 
environments, m geometncally difficult to reach locations, and do all of these at 
relatively large distances from the test object surface Furthermore, lasers are able to 
produce simultaneously shear and longitudinal bulk wave modes, as well as Rayleigh and 
Lamb wave modes LBU system uses a pulsed laser to generate ultrasound and a 
continuous wave (CW) laser interferometer to detect the ultrasound at the point of 
interrogation to perform ultrasonic inspection Some of the important areas of application 
of LBU are 

• Medical field 

• Microwave electronics 

• Defect identification in a wide range of matenals, including metals concrete, 
ceramics and composites etc 

• Aerospace industry 

• Steel industry 

• Nuclear reactors 

• Microstructure analysis of ferrous and non-ferrous matenals 

An additional benefit of laser techniques is that they are capable of a high degree of 
absolute accuracy, since measurements can in pnnciple be calibrated against the 



wavelength of light There may be sensitivity penalties for using lasers rather than 
piezoelectric transducers or probes, and certainly they are likely to be more costly and 
complex to use Nevertheless they are beginning to make a small but significant impact m 
a limited number of applications where their benefits over other probes outweigh their 
disadvantages 

Histoncally, investigators have used the terminology A-, B- or C-scan to descnbe the 
type of data taken during ultrasonic study of specimens An A-scan indicates a variation 
of signal amplitude with time, where the ordinate of the oscilloscope trace is proportional 
to the amplitude and abscissa to the time elapsed A B-scan is the presentation of data in 
the form of a depth profile versus position along a specimen A C-scan is the presentation 
of variation of any feature or feature set as a function of position of transducer in two 
dimensions relative to the surface of the examined specimen Essentially, the feature or 
feature set, used to construct C-scan image is extracted from the received ultrasonic 
waveform obtained m the A-scan mode Another recent way of representing ultrasonic 
data IS by using tomographic technique which reconstructs a slice of an object firom 
multiple projections Thus, tomography, if perfected, can provide better results 


1.2 LITERATURE SURVEY 

Non-destructive evaluation of advanced composite materials poses a challenge to 
researchers and applied technologists Ultrasonic methods, however, aie most common 
and practiced widely because of the low cost, relative ease in their use and the amount of 
information that can be obtained from them In this section, an overall review of 
processing of ultrasomc signals to enhance the confidence of identification of defects m 
composite matenals is presented 

In ultrasonic techniques, information on defect charactenzation requires more evolved 
techniques than classical methods Modem non-destmctive testing of matenals has to 
provide the highest possible detection probability, the correct size and exact onentation 
of defects in the specimen The simplest acoustical charactenzations of matenals deal 



with velocity measurements Some well-known techniques of calculating wave speed by 
measunng the time-of-flight are descnbed in Krautkramer and Rrautkramer [1]. Neslroth 
et al [2] introduced the idea of using vanous features of ultrasonic signal to detect vanous 
types of defects, flaws or anomalies m composite matenals Rose [3] introduced a 
concept of feature mapping The basic concept is based on the fact that different types of 
defect are expected to interact distinctively with different features of ultrasomc signal He 
proposed the Fischer linear discnminant function as a classifier, which uses linear 
combination of useful features to detect anomalies Sann [4] and Dahraju [5] have 
developed an expenmental set-up using an Nd YAG pulsed laser for generation of 
ultrasound m specimens and a He-Ne continuous laser based heterodyne optical 
mterferometnc probe for detection of the same They have proposed a methodology for 
inspection of composite specimens using LBU and evaluation of effectiveness of 
different features of ultrasonic signals in detecting different types of defect Scruby and 
Dram [6] m their introductory book on this subject have elaborated on the pnnciples 
underlying generation and reception processes of all techniques used 

Tomographic imaging pnnciples have been successfully applied to the problem of 
recovenng local matenal properties from acoustic wave propagation data for many years 
While acoustic tomography has many features in common with its radiographic 
counterpart, there are additional complicating features (notably ray bending), which make 
tomography more difficult with acoustic data When anisotropy is included, the level of 
difficulty increases markedly Kline [7] used tomographic imaging pnnciple to 
charactenze the mechamcal properties of anisotropic media His work is mainly focused 
on geophysical applications where weak anisotropy is assumed 

Kline and Wang [8] developed a technique for ultrasomc tomography in anisotropic 
media The ART algonthm was modified, through the use of the governing equations for 
wave propagation in anisotropic media and composite micromechanics Variation of 
acoustic velocity with propagation direction is also taken into account Elline, Wang, 
Mignoga, and Delsanto [9] proposed an approach to ultrasomc tomography in anisotropic 
media based on finite difference formulation Chow, Hutchins, and Mottram [10] 
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conducted expenments in which acoustic emission and ultrasonic tomography were 
performed concurrently on samples of anisotropic glass fiber reinforced polymer 
composite in a direct tensile test Sullivan, Kline, and Mignona [11] developed an 
algonthm for acoustic tomographic reconstraction of fiber angle in orthotropic 
composites The fiber angle about the y-axis was constructed using a dimensional 
solution Rathore, Munshi, and Kishore [12] presented a new tomographic reconstruction 
method, C-ART, for anisotropic matenals, which accounts for the directional dependence 
of sound velocity in the matenal Subbarao, Munshi and Muralidhar [13] studied the 
performance of ART, MART and MENT family of algonthms under limited data 
conditions Their study concludes that MARTS performs better than all other algonthms 
Farhan [14] made a bnef study on ultrasonic tomography of Perspex/polystyrene 
composites, where the different acoustic impedances zones are clearly shown by using 
Convolution Back Projection algonthm 

Ultrasonic measurements are widely used for the elastic constant determination of 
anisotropic matenals Choosing an appropnate propagation direction and polanzation, 
can maximize the accuracy of elastic constant measurement Six elastic constants for an 
orthotropic matenal can be determined from measurements along S 5 mmetry axes 
Degtyar and Rokhlin [15] proposed an inversion procedure based on nonlinear least 
square method to determine elastic constants from group or phase velocity data m 
orthotropic and transversely isotropic matenals Bernard, Jenkins, and Wolfgang [16] 
developed a method, which optimally determines the elastic constants of anisotropic 
solids from wave-speed-measurements m arbitrary nonpnncipal planes Reverdy and 
Audom [17] descnbes a method for determining elastic constants of anisotropic matenals 
from phase velocities of acoustic waves generated and detected by lasers Descharaps and 
Bescond [18] proposed a numencal method to recover the elastic constants from 
ultrasound group velocities 

In the literature, only a few ray-tracing algonthm for wave propagation in anisotropic 
inhomogeneous matenals are discussed Spies [19] studied elastic wave propagation in 
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homogeneous and layered transversely isotropic media using plane wave solutions and 
Gaussian wave packets Rudzki [20] proposed Fermat’s principle for anisotropic media 

1.3 PRESENT WORK 

In the present investigation , time-of-flight (TOF) tomography has been used to 
reconstruct slowness in a glass-epoxy composite domain containing an insert of a 
sphencal steel ball We have taken cross ply composite with insert at center, so symmetry 
in TOF data is expected, and we use this as an assumption to complete the data set ART, 
that IS explained in chapter 3, has been employed for reconstruction which essentially 
computes the slowness distnbution m the domain from a numerous projection TOF data 
from each different view Laser Based Ultrasomcs(LBU) , as discussed above, is used to 
generate A-scan images for excitations along 0, 90, 45, and 135 degrees with respect to 
x-axis TOF information is then extracted from those images One of the objectives of this 
work IS to observe the variation in the response at the receiver end when the concerned 
ray passes through or bypasses the insert region 

The present study has been conducted on a domain (Matenal glass/epoxy) with 
dimensions SOmmxSOmm contammg an inclusion of a sphencal steel ball of 4 76mm 
diameter as shown in fig 4 1(b) For reconstruction purpose, the domain is divided in to 
30x30 gnds and the slowness for each gnd is evaluated in the reconstruction process 
This analysis is repeated for four directions (views), namely 0, 90, 45, and 135 degrees 
with respect to the positive x-axis as shown m the figs 4 l(c)-4 1(f) We have taken four 
A-scan images at 0 and 90 degrees with 3mm ray spacing, while nine A-scan images with 
2 12 mm ray spacing are considered for 45 and 135 degrees Finally, for reconstruction 
purpose 30 rays have been considered for 0 and 90 degree view, while 60 rays are 
considered for 45 and 135 degrees For each view number of rays are chosen such that 
each boundary pixel acts as a source with opposite one as receiver Rays other than those 
from A-scan image are assumed on the basis of symmetry considerations and in-between 
rays are fixed by linear interpolation 
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Numerical Simulation of a similar object is perfoimed to check the accuiacy of the 
various algorithms and to compare with the experimental lesults 


1.4 THESIS ORGANISATION 

The various chapters m the thesis deal with the following aspects 

• Chapter 2 presents the basics of laser ultrasonics 

• Chapter 3 deals with the experimental set-up and procedure 

• Chapter 4 deals with various ultiasonic tomographic reconstruction techniques 

• Chapter 5 presents results and discussion of the present work 

• Chapter 6 finally concludes the present work and gives suggestions for future 
work 
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Chapter 2 

BASICS OF LASER ULTRASONICS 


2.1 INTRODUCTION 


Laser (Light Amplification by Stimulated Emission of Radiation) is a device that 
generates or amplifies light, just as transistors generate and amplify electronic signals at 
audio, radio or microwave firequencies Here, light must be understood broadly, since 
lasers have covered radiation at wavelengths ranging fi’om i n frared range to ultraviolet 
and even soft x-ray range The term ultrasonics is used to descnbe mechanical wave 
propagated m gases, liquids or solids at frequencies above the upper limi t of human 
heanng i e 20KHz and up to 400MHz Because the charactenstics of these waves are 
influenced by the mechanical properties of any medium through which they pass, one can 
use ultrasonics to investigate the properties of that medium 

The laser-based ultrasonics (LBU) is the technique, in which laser beam interaction with 
the surface of the test sample is substituted for piezoelectnc transducers for launching 
and probing elastic waves Since they do not require any mechamcal contact, these 
techniques are very attractive For example, m the field of NDT, the association of laser 
generation with optical detection provides a completely remote ultrasonic system Other 
fields of applications are matenal evaluation, acoustic emission photo thermal 
microscopy and acoustic field imaging According to the intensity of the laser, the impact 
generation method may be classed in two mam categones thermoelastic regime and 
ablation regime fri the field of NDT, surface damage is avoided by means of specific 
techniques In concrete testing, the structure can tolerate some surface imperfection 
There are different methods applied, for example, one uses shock waves, which are 
generated by laser impact on the sample and optical detection (interferometncs or non- 
interferometncs), which monitors the induced surface Other methods are used 
incorporating contact ultrasonic probes working as receivers or transmitters Hence, this 
method needs no contact medium and the source can be some meters away, it is used 



especially for applications such as NDT of systems at high temperature 

2.2 MAIN PROPERTIES OF LASER LIGHT 

Lasers are charactenzed by a number of key optical properties, most of which play an 
important role in the interaction with ultrasonic fields Their four major optical properties 
are 

1 Monochromaticity 

2 Coherence 

3 Directionality 

4 High Intensity or Brightness 

2.2.1 Monochromaticity 

Monochromatic means the same frequency Light from sources other than lasers covers a 
range of frequencies True single-frequency operation of a gas laser can be achieved by 
careful design Monochromaticity is reduced by multimode operation of a laser A small 
He-Ne laser generally has three or four longitudinal modes excited with a spacing of a 
few hundred MHz, depending on the cavity length This still gives a bandwidth of a few 
parts m 1 0 ^ Solid-state lasers tend to have rather large frequency spreads 

Monochromaticity is important for some ultrasonic applications, m particular the 
mterferometnc measurement of ultrasonic fields Firstly, it is necessary in order to obtain 
high coherence, and secondly, it enables accurate absolute measurements of ultrasonic 
displacement to be made, since calibration is against the wavelength of the light 
Monochromaticity is of less importance for the generation of ultrasound by laser 

2.2.2 Coherence 

Coherence means the same order or as a copy of the other photon and in phase with the 
other photon Coherence is an important property when it comes to building an optical 



system to detect ultrasonic waves In simple words, coherence is used to descnbe how 
well a wave disturbance at one point in space or time correlates with the disturbance at 
another point If there is a well defined phase relationship between the light at two 
different points in space (i e two light beams), or at two different times (i e one beam 
split into two with a delay between the parts, as is typical in an interferometer) then the 
two light disturbances can be brought together to produce a predictable interference 
pattern If the light is completely incoherent so that there is no predictable relationship 
(i e random phase) between the two disturbances, no interference fhnges will be formed 
Fig 2 1 shows the comparison between the ordinary light and laser hght 

Coherence is important to the reception of ultrasonic waves by laser Most techniques 
involve some form of interferometer, in which good coherence between probe and 
reference beam is essential Coherence length is the distance of the ongin of the beam to 
the farthest point at which the wave disturbance can be effectively correlated with the 
disturbance at its starting point Conventional monochromatic light cannot be used 
because their coherence length is only of the order of millimeters The use of a gas laser 
however enables use of a longer probe than reference beam This means that a compact 
instrument (incorporating all except the probe beam) can be built, and that the distance to 
the sample is not cntical 

2.2.3 Directionality 

Directionality is a function of the spatial coherence of the beam of hght Laser beam is 
highly directional which implies laser light is of very small divergence because of the 
parallelism of the beam The radiation produced by a laser is confined to a narrow cone 
of angles The beam of divergence for a typical gas laser is of the order of 1 milhradian, 
although some commercial He-Ne lasers are available with divergences of a few tenths of 
a milhradian Comparable beam divergences are now available fi'om pulsed solid-state 
laser systems It is only possible to attain such low beam divergence from conventional 
hght sources by severe colhmation, which leads to extremely low intensities The 



diameter of the beam is typically about 1mm for a gas laser such as helium-neon, and in 
the range 1-20 mm for a pulsed solid-state laser 

The light from gas and solid state lasers thus forms a highly collimated beam which is 
extremely valuable for laser ultrasonics since it enables the beam to be focused to a very 
small spot This means not only high spatial resolution, but also, in the case of ultrasonic 
displacement measurement by interferometer, the ability to collect a larger fraction of the 
scattered light from a rough surface, thereby increasing sensitivity For laser generation it 
indicates that very high incident power densities are attainable Low beam divergence 
also means that the beam can travel distances of the order of several meters from the laser 
to the specimen without appreciable spreading and losses Thus both laser generation and 
reception of ultrasound can be made genuinely remote techniques. 

2.2.4 High Intensity or Brightness 

Intensity or bnghtness of a light source is defined as the power emitted per unit area per 
unit solid angle A laser beam is extremely intense, more intense than any other light 
source This is perhaps the property for which lasers are best known outside the field of 
optics Although the optical power output from a small helium-neon (He-Ne) laser may 
only be say 2m W, a beam diameter of 0 5mm leads to a power density of about 1 W cm'^ 
Such a beam can readily be focused by a simple lens to a spot of diameter 0 05 mm 
because it is monochromatic and coherent The incident power intensity is then 100 W 
cm"^ Focusing optics produces a beam of sufficient intensity to melt, cut or weld 
structural matenals 

Intensity is a most important property for laser reception of ultrasound, as sensitivity of a 
single mode laser mterferometer system (defined as signal to noise ratio for a fixed 
bandwidth) increases with the square root of light intensity, provided other conditions 
remain constant The limiting factor becomes the intensity at which the specimen is 
damaged or otherwise adversely affected by intense irradiation Increase in power also 
bnngs penalties like increase m noise and multimode interference Intensity is also a 



crucial factor in the generation of ultrasound by laser since incident power intensities 
typically in the range of 10"^- 10^ Wcm'^ are needed to act as a theimoelastic source of 
ultrasound The laser energy per pulse largely controls the ultrasonic amplitude in the 
thermoelastic regime The limiting factor is not the power that is available from 
commercial lasers, but rather the threshold for damage m the irradiated sample. Indeed, 
many pulsed lasers are too powerful for ultrasonic applications, and must be used at 
reduced power if the technique is not to be destmctive 

2.3 BASIC PRINCIPLE OF LASER OPERATION 

The laser is a device, which amplifies the intensity of light by means of a quantum 
process known as stimulated emission A practical laser needs a means of amplification 
(i e stimulated emission) and a means for feeding the energy back into the system to 
build up sustained oscillation The operation of the simplest type of laser can be most 
readily understood in terms of a quantum mechamcal model having say three energy 
levels Eo, Ei, E 2 , such that E 2 >Ei>Eo (Fig 2 2) In reality, there may be more than three 
levels, but for present purposes, other levels are neglected The ground state Eo is well 
populated, whereas the intermediate and upper states are more sparsely populated 
Suppose now the atom absorbs a quantum of incident radiation such that the atom is 
excited into the upper state From quantum theory, the radiation must have a frequency Vp 
such that 

hVp = E2-Eo (2 1) 

Where h is the Planck’s constant In laser terminology this process of absorption is 
known as ‘pumping’, so that Vp is the ‘pumping frequency’ Pumping tends to equalize 
the population of two states so that E 2 becomes well populated This reverses the normal 
occupancy of E 2 and Ei and is known as population inversion Emission (1 e stimulated 
emission) can now occur in response to incident radiation, at a frequency v given by 

hv = E 2 - El (2 2) 

It can be noted that necessanly Vp > v, so that the pumping frequency must always be 
equal to or higher than that of the radiation to be amplified In order to obtain light of 



sufficient intensity for practical use, there has to be some mechanism for feeding the 
energy back into the laser system and thereby building up the amplitude of oscillations m 
a resonant system The usual way of obtaining sustained oscillations is to site a high 
performance mirror at each end of the lasing medium In the simplest system, both 
mirrors are plane, and accurately aligned perpendicular to the axis of the laser Thus, the 
light can be reflected backwards and forwards through the lasing medium On each pass 
it stimulates further emission firom the medium and is thus amplified in intensity 

2.4 LASERS FOR ULTRASOUND GENERATION 

The simple laser arrangement will operate m a pulsed mode if it is pumped for example 
by a pulsed flash tube Depending on the type of laser, pulses of duration typically 100 ps 
to 1ms can be obtained Although high-energy pulses can be produced m this way, the 
normal mode is not particularly useful for laser ultrasonics because the pulse duration is 
too large An additional technique, known as Q- switching or Q-spoiling, is needed to 
obtain pulses m the required 1-1 00ns range The Q (quality) factor of a cavity resonator is 
the energy stored m the cavity divided by the energy lost from the cavity per round tnp of 
the light within the cavity Thus, if Q is low, the cavity oscillations are suppressed and 
the stored energy builds up within the lasing medium When the Q is high, the cavity can 
support oscillations mto which energy is supplied from the medium Thus switching from 
low to high Q results m the rapid extraction of power from the laser cavity Practical Q- 
switches are m the form of elements with variable absorption that are inserted between 
the mirrors Two commonly used Q switches are the Pockels cell and bleachable 
(saturable) dye Also of potential importance to non-contact ultrasonics are lasers that can 
be pulsed repeatedly Adequate cooling m pulsed solid-state lasers is a problem Hence 
gas lasers are preferred However solid-state lasers can store higher mdividual pulse 
energies than gas medium lasers Conventional ultrasonic inspection often uses repetition 
rates as high as 1-10 kHz for signal averaging purposes or for speed 

Laser generation of ultrasound does not impose very stnngent conditions on the source 
laser All that is required is the ability to deliver a reasonably high-pulsed energy density 



to a small area of the specimen, where the pulse length is in the range of 1-100 ns As 
already mentioned, wavelength is not cntical, nor is monochromaticity or coherence 
Most researchers have employed a solid state laser, either mby or Nd YAG The Nd 
YAG laser has proved itself to be a versatile system under Q-switched conditions 
Adequate pulse energy can be obtained from Nd YAG lasers of modest size, while the 
pulse lengths are close to ideal for ultrasonic studies m metals However, the fundamental 
wavelength (1064nm) is in the near infrared, which makes alignment more difficult than 
with a visible laser, in addition to increasing the nsk of accidental eye damage The Nd 
YAG laser is available with frequency doubling to give several harmonics of the 
fundamental frequency The frequency-doubled wavelength of 532 nm is m the middle of 
the visible spectrum, while the harmomcs of 266nm and 355nm are in the ultraviolet 
zone 


2.5 LASERS FOR INTERFEROMETRY 

For applications like interferometry, where wavelength punty and coherence are 
important, continuous wave lasers are to be designed with special optical components 
The mam problem is that a simple laser will excite a number of longitudinal and 
transverse modes The result is that energy is amplified over a narrow range of 
frequencies instead of the desired single frequency, giving nse to intermode beat 
frequencies, and there is a vanable distribution of energy across the beam The generally 
preferred transverse mode (usually the lowest order mode with circular symmetry, 
TEMoo) can be selected, and higher order modes suppressed, by the use of at least one 
curved (concave) mirror at the end of the cavity and/or a suitable aperture within the 
cavity A multi-mode laser can be employed to mimmize the multi-mode problem of 
longitudinal wavelengths in interferometry and then a balanced detection system and 
adjustment of the path difference between the two arms of the interferometer can be used 
to minimize the effects of intermode beats if necessary 

The reception of ultrasound by laser is dependent upon all the fundamental properties of 
the laser monochromaticity, coherence, directionality and high power density Because 



of its ready availability and excellent optical charactenstics, notably monochromaticity 
and coherence, the He-Ne laser has dominated much of the work This system is however 
somewhat limited with regard to maximum power which controls sensitivity The other 
choice IS the argon ion laser, which not only can deliver higher power, but also inherently 
more sensitive because of its shorter wavelength Argon ion lasers are however very 
bulky and noisy The CW Nd YAG laser is also available as a solid-state option, 
although it suffers from a number of drawbacks including bulk 

2.6 LASER-IRRADIATED ULTRASOUND MODES 

Laser ultrasomcs are based on the phenomenon that laser irradiation onto a solid surface 
generates ultrasonic waves m the solid Fully non-contact generation and detection of 
ultrasonic waves can be realized by combimng with optical (laser) interferometer for 
ultrasomc detection Thus, laser ultrasomc NDE could be applied to matenals and 
structures m hostile environments, such as very high temperatures, or to components 
having complex geometnes A laser emits a beam of coherent radiation, whose 
wavelength may be in the infrared, visible or ultraviolet part of the electromagnetic 
spectmm When this is incident on a solid sample, in general some of the energy is 
absorbed by various mechamsms, depending upon the nature of the sample and the 
frequency of the radiation, while the remainder is reflected or scattered from the surface 
Three basic mechanisms have been illustrated below 

(I) The Thermoelastic Regime: At lower incident powers, electromagnetic radiation 
from the laser is absorbed in the surface region of a sample, causing localized heating as 
shown m Fig 2 3 Thermal energy then propagates into the specimen as thermal waves 
The heated region undergoes thermal expansion and thermoelastic stresses generate 
elastic waves (ultrasound), which propagate deep within the sample For typical Q- 
switched laser pulse durations, the thermal wave field only extends a few micrometers 
even m good conductors Contrast the incidence of low frequency modulated light where 
the thermal field extends millimeters or centimeters and is itself useful for matenals 
charactenzation 



(II) The Constrained Surface Source- Constrained surfaces include coating the surface 
with a thin solid layer of, for instance, paint or rust and roughness, covenng the surface 
with a transparent solid such as glass, covenng the surface with a transparent liquid, and 
finally constraining a thin layer of liquid between a transparent solid and the sample All 
these four techniques for modifying the surface of the sample introduce large stresses 
normal to the surface, as shown m Figs 2 4 (a-d), which are otherwise absent for the 
thermo elastic source at a free source The form of the ultrasomc source is thus 
substantially modified The different combination of stresses generates a different 
ultrasonic sound field, with considerable enhancement m compression-wave amplitude, 
especially propagating normal to the surface The constrained surface source generates a 
sound field more akin to that of ablation, but at lower power densities and without any 
damage However, the need to constrain with for instance a liquid, reduce the usefulness 
of the laser ultrasonic source for many non-destructive applications 

(III) The Plasma Regime: At higher incident powers, surface melting and evaporation 
occurs, resulting m material ablation and the formation of plasma above the sample 
surface as shown m Fig 2 5 The momentum of the evaporated matenal exerts an opposite 
force on the sample, causing a reactive stress at the surface This generates an intense 
broadband ultrasomc source, whose strongest components are, directed normal to the 
surface The effects caused by the plasma generated in association with ablation in case 
of Q-switched lasers are 

• The plasma exerts a high pressure on the surface, which in turn suppresses 
vaponzation of the matenal by raising the boiling point of the matenal well above 
Its normal value 

• It absorbs light from the laser pulse, acting as a shield, but also becoming 
extremely hot 

• As it expands it produces an impulse reaction on the surface 

• It radiates some of its heat back on to the surface, maintaining its high 
temperature for some time after the incident laser pulse power has started to fall 




Fig 2 1 Companson of ordinary light and laser light 
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Fig 2 3 Thermoelastic regime 
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Fig 2.4 Schematic diagram showing stresses induced when a laser pulse is incident on a 
sample surface covered by (a) a thin layer of paint, (b) a thick transparent solid (c) a layer 
of liquid (d) a layer of liquid constrained by a transparent solid 
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Fig 2 5 Schematic diagram to show ablation of surface matenal and net reactive force on 
sample 

2.7 PRINCIPLE OF LASER INTERFEROMETRY 

Optical interferometry is a very sensitive way of measunng displacement, but to be 
practicable for general use, it requires a highly monochromatic light source and thus the 
use of lasers is virtually essential Interferometers for the detection of ultrasonic 
movements of waves may be divided into two mam types In the first type, light scattered 
or reflected from a surface is made to interfere with a reference beam, thus giving a 
measure of optical phase and hence instantaneous surface displacement The second type 



of interferometer makes use of interference between a large number of reflected beams 
This IS designed as a high-resolution optical spectrometer to detect changes m the 
frequency of the scattered or reflected light It thus gives an output dependent on the 
velocity of the surface The first type is the more widely used and the most practical at 
lower frequencies and with reflecting surfaces The second type offers a potentially 
higher sensitivity with rough surfaces, particularly at higher frequencies 

For the detection of ultrasomc waves at a surface, the techmques are admittedly 
insensitive compared with piezoelectnc devices They do, however, offer a number of 
advantages 

(1) They are non-contacting and thus do not disturb the ultrasomc field The point of 
measurement may be quickly moved and there are no fundamental restnctions 
on surface temperature 

(2) High spatial resolution may be obtained without reducing sensitivity The 
measurements may be localized over a few micrometers if necessary 

(3) As the measurements may be directly related to the wavelength of the light, no 
other calibration is required 

(4) They can have a flat broadband frequency response, something difficult to 
achieve with piezoelectnc transducers, particularly at high frequencies 

2.8 HE-NE LASER HETERODYNE INTERFEROMETER 

2.8.1 Working of Optical Heterodyne Probe 

The system consists of a compact optical head and an electronic signal-processing umt 
The optical layout is shown in Fig 3 2 The laser beam emitted by the laser source, 
honzontally polarized, is directed in the interferometer by two deflecting mirrors The 
laser beam is split into two parts by the beam-splitting cube The reference beam is 
reflected by the beam-splitting cube and goes through the Dove pnsm It is then 
transmitted by the polanzmg beam-splitting cube and deflected by the mirror on the 
photodetector The probe beam is transmitted by the beam-splitting cube and is 



horizontally polanzed Its optical frequency is shifted in the Bragg cell, and transmitted 
by the polanzmg beam-splitting cube A quarter wave plate transforms the honzontal 
polanzation into circular polanzation The lens focuses the beam on the surface of the 
sample The probe beam is then phase modulated upon reflection on the sample by the 
mechanical displacement After the second pass in the quarter wave plate, the direction of 
polanzation becomes vertical The probe beam is reflected by the polanzmg beam-splitter 
on the photodetector Just after the polanzmg beam-splitter, the polanzation of the probe 
and reference beam are respectively vertical and honzontal These two beams can 
therefore not interfere The analyzer selects a common component at 45° of the two 
polanzations, thus allowing interference The photodetector delivers a beat signal at the 
frequency of the Bragg cell phase modulated by the mechanical displacement of the 
object 


2.8.2 Principle of Detection 

The complex amplitude of a laser beam of frequency frcan be wntten as 

L = (2 3) 

This is divided in the interferometer into a reference beam and a signal beam, whose 
complex amplitude is 

i? = r (2 4) 

The reference beam does not expenence any perturbation, however the signal beam 

expenences a frequency shift fs m the Bragg cell Upon reflection on the object, its phase 

is modulated by the displacement of the sample 

(p(t) =471 d(t)/A, (2 5) 

where, X is the wavelength of the laser beam and d (t) the mechanical displacement of the 

object The complex amplitude of the signal beam is therefore 

^ ^ g2in/ii+2iri/si+/p(0 ^2 5) 

The interference of the two beams on the photodetector produces an electncal signal at 

frequency fe, phase modulated by the displacement of the object 

I(t) = Io + i(t) (27) 



I (t) = kcos (27i:fBt + {p(t)) (2 8) 

The useful signal is contained in the signal delivered by the photodetector as a phase 
modulation of the earner frequency The signal processor delivers an electnc signal 
proportional to the displacement of the object Half of the current i(t) is filtered at the 
frequency fg, and phase shifted by 90° It is then mixed with the other half, non 
perturbed, and yields a current 

j (t) oc cos (27tfBt + (p(t)) X cos (27i:fBt + 7t/2) (2 9) 

j (t) oc l/2[cos (47rfBt + 9 (t)+ Tt/2)+ cos ((p(t)+ 7r/2)] (2 1 0) 

The signal at the frequency 2fB is filtered, to give 
s (t) oc sin (p(t) (2 1 1) 

If the displacement is very small compared to the optical wavelength, this signal can be 
wntten as 

s (t) = k 471 d (tyX (2 12) 

The final electneal signal is therefore directly proportional to the displacement of the 
object 


2.8.3 Sensitivity 

Laser ultrasonics suffers from a lack of sensitivity relative to conventional ultrasonics as 
the photon structure imposes a fundamental limit on the change m light levels dunng 
generation / detection of lasers 
Overall sensitivity to flaws is given by 

Sensitivity=T x /(ct. A) x R (2 13) 

Where T, /(a, A), R are terms related to the transmitted sound, the focal or imaging 
properties of the system determined by the flaw scattering a and the focal aperture A, and 
the receiver sensitivity The generated signal can be improved using tailored surface 
coatings or by increasing the power of the generating source However, this approach is 
limited by the onset of ablation m the target 



2.9 CLOSURE 


In this Chapter, the basic theory of geneiation of ultrasonic waves by lasers and their 
detection using laser-based interferometers has been presented Some of the important 
properties of lasers related to aspects of generation and detection of ultrasonic waves 
have been discussed Advantages of using laser-based detectors despite poor detectivity 
and sensitivity have also been explained 



Chapter 3 

EXPERIMENTAL SET-UP AND PROCEDURE 


3.1 INTRODUCTION 

This chapter bnefly descnbes the details of laser-based ultrasonic (LBU) set-up, 
developed for non-destructive evaluation (NDE) of composite matenals and data 
acquisition procedure 

Nd YAG pulsed laser is used to generate ultrasonic waves in matenals or objects to be 
inspected Heterodyne type laser interferometer is used to detect the transmitted wave 
through the matenal, the signals are then amplified and digitized using a Yokogawa 
DL1740 Digital Storage Oscilloscope (DSO) The oscilloscope is tnggered using a 
synchromzation signal from the pockels cell of the pulsed laser The recorded waveforms 
through DSO are transferred to a HP Pentium IE computer over an USB/Ethemet 
interface for subsequent storage and analysis The schematic layout of expenmental set- 
up IS shown in Fig 3 1 and photograph of the same is shown in Fig 3 3, while Optical 
Layout of Heterodyne Interferometer is shown m Fig 3 2 The scarming is done manually 
using two single axis micrometer controlled XYZ translator mounted on the Optical Test 
Bench A bnef descnption of the components of the present set-up is given in the 
following paragraphs 

3.2 EXPERIMENTAL SET-UP: LASER-BASED ULTRASONICS 

3.2.1 Nd; YAG Pulsed Laser Ultrasonic Generator 

The 5000 DNS senes pulsed Nd YAG laser is built on a modular concept The optical 
configuration allows vanable setting of the optical parameters It consists of three major 
components, the laser-head, the power supply and the cooling unit The heart of the 



system is the pumping stmcture, which houses the Nd YAG rod and the flash lamp The 
lasers are built on an elect o-optically Q-switched oscillator This oscillator uses a 
pockels cell Q-switch to produce pulses of high intensity and short duration (i e m the 
order of 5-7 ns) Q- Switched lasers are often used as a non-contact ultrasound source m 
non-destructive testing (NDT) of matenals Q-switched lasers typically have nanosecond 




(ns) pulse durations and generate broadband ultrasound waves, though longer laser pulses 
of 100 microseconds or greater have been used for NDE A vanable reflectivity output 
coupler allows the extraction of high energy on a single spatial transverse mode This 
leads to laser beams of low divergence, and to high conversion efficiencies in the 
harmonic wavelengths (1064, 532, 355, 266 run) Different harmonic generators extend 
the wavelength range to the second, third and fourth harmonics 

The active medium of the laser operates on transitions of tnply ionized Neodymium 
atoms (Nd3+), which take place of another ion (Yttnum) in the host Yttnum Aluminum 
Garnet known by the acronym YAG The laser operates as a 4-level system The intense 
broad-spectrum light of the flash lamp populates the upper level Once m the higher 
energy level, the Neodymium ions drop to a metastable level, producing a population 
inversion The lower level decays by a fast non-radiative process to the ground state The 
strongest Neodymium line is 1064 nm 

3.2.2 Optical Heterodyne Laser (He-Ne) Probe 

The SH-130 probe is designed to measure transient mechanical displacements of very 
low amplitude It is specially devoted to measure displacements generated by the 
propagation of an acoustic or ultrasonic wave The system consists of a compact optical 
head and an electronic signal-processing unit The optical head integrates high stability, 
low power laser source for fast detection with a high spatial resolution 

The electronic signal processor delivers a response proportional to the displacement of 
the target, with a high bandwidth The output signal is automatically calibrated to give the 
absolute value of the measured displacement The pnnciple of detection (heterodyne 
interferometry) makes the system insensitive to external vibrations The compactness of 
the system allows for a wide range of operating conditions 




Figure 3 2 Optical Layout of Heterodyne Interferometer 


3.2.3 Digital Storage Oscilloscope 

The set-up utilizes a Yokogawa DL1740 (four channel, one GSa/sec, 500 MHz) Digital 
Storage Oscilloscope, having bmlt in Zip dnve, Ethernet, USB, GPIB and Serial Ports for 
communication wath external PC’s/ systems 










Figure 3 3 Photograph of expenmental set-up 

3 2 4 Computer 

The DSO IS interfaced to a Pentium III based PC and the data is transformed on-line 
using communication ports (Ethernet) 
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3.3 EXPERIMENTAL PROCEDURE (LBU) 


3.3.1 Calibration of Set-up 

The system was checked for calibration before scanning the specimen The procedural 
steps involved to check for calibration are as follows 

(1) Alignment of the SH-130 Optical Probe This involved optimizing the contiol 
signal by adjusting the onentation of highly reflective mirror and fine-adjustment 
of its distance fiom the focusing lens (The focal length of the focusing lens is 
around 215 mm ) 

(2) The photo detector output was measured without going through the signal 
processor This was found to be 310 mV (as per specifications) The DSO image 
IS shown in Appendix ‘A’ 


(3) The output was then taken through the signal processor with the automatic gam 
control switched on The measured signal was found to be in excess of the 
stipulated signal level of 630 mV peak to peak on 50Q (found to be 730 mV) 
This was adjusted using the adjustments available m the signal-processing unit 
The DSO image is shown in Appendix ‘B’ 


3.3.2 Preparation of Specimens 

The present ultrasonic non-destructive evaluation was done on composite specimens 
implanted with artificial inclusions like metal inserts The specimens were prepared by 
hand lay up technique, using woven glass fabric as the reinforcing material and epoxy 
resin as the matrix material Resin used being LY556, hardener HT976 and accelerator 
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XY73 all manufactured by Ciba Geigy (India) The curing was done at 120° C and 5 atm 
pressure for one hour followed by one hour at 150° C and 5 Atm pressure and left to cool 
to room temperature at same pressure overnight 

The dimensions of the finished specimens were 30 x 30mm^ Model specimens were 
prepared with a central mild steel insert to simulate presence of flaws or defects The 
metal inserts, spherical in shape, about 4 76mm m thickness was introduced centrally in 
the thickness direction (m between plies 36 and 37) All the specimens were cast with 72 
plies with a final thickness of around 20 mm 

3.3.3 Scanning Procedure 

The specimen to be investigated was cleaned and the area under investigation was 
marked To get good signals, the portion of the specimen facing to the Ne-He laser was 
pasted with retro-reflective tape The He-Ne laser was switched on to align the CW laser 
beam with the marked point of investigation to obtain good information regarding the 
wave propagation Then the Nd YAG laser was focused onto the specimen using an 
arrangement of optical mirrors and lenses (as shown in figure 3 1) so that it acts as a point 
source and is perpendicular to the irradiated face of the specimen The alignment of Nd 
YAG and He-Ne laser was ensured by focusing the Nd YAG laser on the marked point 
of investigation on the opposite face of the specimen The power level of the laser was 
adjusted to ensure thermoelastic region so that there was no ablation of the specimen The 
scanning was performed at setting corresponding to 1064nm wavelength for the Nd 
YAG laser The He-Ne laser was focused on the opposite side of the specimen at the 
focal length of the collection lens (approx 215mm) The output of the Optical probe was 
fed to the electronic signal-processing unit The electronic signal processor delivers a 
response proportional to the displacement of the target, with a high bandwidth The 
output signal was automatically calibrated to give the absolute value of the measured 
displacement This was fed to the Digital Storage oscilloscope (DSO) where it was 
displayed and storage of the data was carried out on-line m the PC interfaced to the DSO 
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3.4 CLOSURE 


In this Chapter, the details of experimental set-up and scanning procedure are presented 
for the Laser-Based Ultrasonic (LBU) set-up Experimental procedure for data 
acquisition is also explained 
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Chapter 4 

ULTRASONIC TOMOGRAPHIC 
RECONSTRUCTION 


4.1 INTRODUCTION 

The technique of Computenzed Tomography (CT) has established itself as a leading tool 
in diagnostic radiology over the past twenty years and it is gaming importance m the 
nondestructive evaluation area m a vanety of situations In conventional projection 
radiography, the result is just obtained just in the form of a shadowgraph reqmnng a lot 
of interpretation CT, however, uses tomographic reconstruction to combine the 
information from multiple projections Thus, the CT information is displayed 
quantitatively as a reconstracted slice plane of the object and it can provide accurate 
quantity measures of matenal charactenstics in small volume elements of the component 
CT was ongmally conceived for X-rays, but m recent times, CT has been implemented 
for other energy sources such as acoustic, laser, microwaves etc Among these, acoustic 
tomography (especially using ultrasonic waves (UCT)) has special importance m 
application to NDE Like X-ray CT, the physics of ultrasomc CT involves two steps 

1 The measurement of the time of flight data for the wave passing trough a number 
of different paths through a test object, followed by 

2 The reconstruction of slowness coefficient values at all the interior points from 
the data collected m step (1) 

Thus, acoustic tomography can be defined as a technique to determine values of a 
spatially varying acoustic parameter of interest across a slice through an object Waves 
travel through the sample (from a senes of source locations) to receivers at appropnate 
positions around the object The source and the receiver locations are chosen such that 
the rays pass through a large portion of the object volume as possible The relative 



positioning of the source and receiver should also conform to any requirement for regular 
positioning m the reconstruction procedure As the source and the detector are moved 
around the object, a number of individual projections (or measured data sets) are 
obtained Each projection is composed of a transmitted intensity pattern, which depends 
on the path concerned After completmg a scan consisting of a number of projections, the 
distribution of the slowness coefficients withm the object is reconstructed by one of the 
several reconstruction techniques 

The major problem in applying UCT is that ultrasound does not travel stnctly in linear 
fashion and suffers scattering, ray bending and the diffraction at the interface of any 
dissimilar matenal In medical applications, however, UCT can be used safely in soft 
tissues imaging where deviation from linearity in propagation is imnimal Among 
industrial applications, practical matenal charactenzation problem such as study of 
residual stress distnbution has also been tackled by UCT With proper guideline for 
collecting projection data and proper modification in the reconstruction technique, 
accurate quantitative evaluation of the matenal properties of the object slice can be 
accomplished 

4.2 RECONSTRUCTION METHODS 

Tomography can be classified into (a) transform (b) senes expansion and (c) optimization 
methods Transform methods m general, require a large number of projections for a 
meaningful answer However, in most of the cases it is not practicable to record such a 
large number of projections partly owing to inconvenience and partly due to time and 
cost Hence, as a rule requirement of a large number of projections is not desirable and 
one must look for methods that converge within a few iterations Limited-view 
tomography is best accomplished using the senes expansion method As limited data 
tomography does not have a umque solution, the algonthms are expected to be sensitive 
to the imtial guess of the field that starts the iterations Optimization-based algonthms are 
known to be independent of initial guess, but the choice of optimization function plays an 
important role in the result obtained Depending on the mathematical defimtion used, the 



enliopy cxlremizalion route may yield good results, while the entropy minimization 
principle may be suitable m other applications For the algebraic technique considered m 
the present study, an unbiased initial guess such as a constant profile was seen to be good 
enough to predict the correct field variable Tomography, being an inverse technique was 
seen to preserve the noise in the initial data However, the dominant trend in the field 
variable was captured during tomographic inversion Actually the problem of 
reconstruction is a problem of inversion of the rectangular matnx Iterative techniques 
that are used in tomography can be viewed as developing a generalized inverse of the 
matnx The general purpose matnx libranes cannot be used to invert them since they are 
highly ill conditioned and rectangular in stmcture The tomographic algonthms can be 
seen as a systematic route towards a meaningful inversion of the matnx 

4.2.1 Algebraic Reconstruction Technique (ART) 

ART IS one of the most popular techniques based upon a row action method for solving 
large, sparse systems It assumes a solution for the image matnx and iteratively compares 
projections through the image to the measured projection data With each iteration the 
image matnx is altered until a match between the calculated and measures projections are 
within acceptable accuracy This techmque is useful when limited angles of projections 
are available and can become highly computation intensive 

In tomography, projection data (p) represents an integrated information resultmg firom 
interaction of source radiation with the matenal, and is expressed as 
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Figure 4 1 Weight function fori*'’ ray and j*'’ pixel 

The iterative reconstruction methods of tomography are based on the descretization of the 
cross sectional plane by a square gnd as shown m the Fig 4 1 The length of the mtercept 
of an i"' ray with the j"' pixel in a given projection is called the weight function (d^ ) If 

IS the field value in the j"' pixel, the i"' projection is given as 

,M (4 2) 

7=1 

where M is the total number of rays and N is the number of pixels 

It IS required to collect the projection data over a number of view angles to carry out 
tomographic reconstruction When ultrasound wave is used as a source of radiation, f(x, 
y) and are usually taken to be the slowness field (inverse of group velocity) in the 
object and time-of-flight (TOF) respectively The aim in tomography is to obtam f(x, y) 
from given p. The conventional algebraic reconstruction techmque (ART) involves four 

major steps 

1 hntial assumption of the object field, 

2 Calculation of correction, 

3 Application of the correction and 





4 Test for convergence 


Typically, m Gordon ART, the field values in pixels are updated for k*'’ iteration using 
concction term as 


/,*=/,*■'+ 




D.. 


D' 


Here the inner product notation is used in the following sense 


Where 




7=1 


A-1 


and 


D' 




(4 3) 


(4 4) 


(4 5) 


(4 6) 


To tackle the problem for inconsistent system of equations, a modification m the ART 
discussed above, is necessary The behaviour of ART algorithm for inconsistent systems 
when relaxation parameters are allowed was discussed by Eggermont et al [1981] 
Relaxation parameters area sequence of real numbers usually confined in some interval 
The modified step (with relaxation parameter) for updating the iterate can be expressed 





(4 7) 


where 

e^<A^<2-£2 ei,£2>0 


Following are some of the other ART algorithms 
• Simple ART 



It IS similar to Gordon ART, except the modification step (with relaxation 
parameter) is given as 
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. Gilbert ART 

Gilbert (1972) has developed independently a form of ART called the 
simultaneously iterative reconstmction techmque (SERT) hi SIRT, the elements 
of the field function are modified after all the correction values corresponding to 
individual rays have been calculated The algonthm is similar to ART but the 
correction is applied as given below 

* Identify all the rays (A^) passing through a given cell and the corresponding 
D‘ and p, 

* For each cell j apply the algebraic sum of all possible correction terms as 




IC=1 


p. - 


ly 


IcA 


D„ 


(4 9) 


4.2.2 MART ALGORITHMS 

The difference between ART and MART algorithms lies in the manner in which the 
corrections are applied to guessed field values The corrections are additive in former, 
while multiplicative in the latter These algorithms are essentially iterative in nature, but 
intermediate steps may also involve repetitive calculations To identify the iterative loop, 



stai t and close labels with statement numbers have been indicated in the description of 
each algorithm 

The MART algorithms are first summarized Let be the projection due to i^'^ray m 0 

direction of piojection and /, be the initial guess of the field value The projection 
using the current field value can be obtained numencally as 

= E A,./, .Me (4 10) 

/=1 

where Me represents number of rays along the 0 direction The parameter to be varied m 
each loop is indicated m brackets 

Start 1 Start iterations (k) 

Start 2 For each projection angle (0) 

Start 3 For each ray ( ) 

Compute the numencal projection [Eq 4 4] 

Calculate the correction as 

Start 4 For each cell ( j ) 

If Wi 0 j IS nonzero then 

MARTI 

x(10-f,x(l- b,it>J) 

MART2 

fr=ff >‘1 

MARTS 

jjiew_j-old 

where p is relaxation factor 
close 4 


(411) 

(4 12) 

(4 13) 

(4 14) 



close 3 


close 2 

Check for convergence as 
If 

loose (415) 

wheic c IS a suitable stopping ciiterion STOP 
Else Continue 
Close 1 

Steps 3 and 4 form the essence of the reconstruction algorithm All three versions include 
the relaxation factor p Typical values foi relaxation factor reported are in the range 0 01- 
0 1, largei values leading to diveigence It is to be noted that correction calculated in step 
3 is the ratio of the recorded projection data ) and that calculated from guessed field 

namely, , which is being iterated The three versions of MART are In MARTI, the 

weight function is prescribed in binary form, being unity if a pixel ray passes through a 
pixel and zero otherwise In MART2 and MARTS, the weight function is calculated 
piecisely as the latio of the length of the ray intercepted by the pixel and the maximum 
dimension of pixel enclosed m it 


4.3 Reconstruction errors 


Dining the picsent woik, the peiformance of leconstiuction technique is evaluated both 
quantitatively and qualitatively The average ELI and RMS EL2 error are used to 
measure performance and defined as 
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4.4 CLOSURE 


In this chapter vanous tomographic algonthms were discussed Error estimation of the 
independent parameters to evaluate the performance of the algonthm was included 



Chapter 5 
RESULT AND DISCUSSION 


5.1 INTRODUCTION 

Firstly, analysis of the signals (A-scan) obtained from Laser Based Ultrasonics (LBU) 
set-up IS done to extract the Time of Flight (TOF) information for subsequent 
reconstructions Then, the leconstruction of slowness field of the crossection under 
consideration is performed using both ART and MART for different values of relaxation 
parameters and initial field values 

As mentioned earlier, a Cartesian grid of square picture elements, called pixels, is 
introduced into the domain (30mmx30mm) of interest Number of rows and columns m 
the grid are 30 Thus the domain is essentially discretized into 900 pixels The numbering 
airangements of the pixels, columns and rows have been chosen m the manner as shown 
m fig 5 1(a) The base material of the domain is glass/epoxy and it contains a cross 
section of steel ball as shown m fig 5 1(b) 

The reconstruction is conducted from projection data Projections have been taken in four 
diiections oi views Four directions, namely 0 90*’, 45°, and 135° with the x-axis, of the 

views or rays have been used m this work 

Finally, for reconstruction purpose 30 rays have been considered for 0 and 90 degree 
view, while 60 rays are considered for 45 and 135 degrees For each view number of rays 
are chosen such that each boundary pixel acts as a source with opposite one as receiver 
Source is considered to be the center of the boundary side of the pixel, so ray spacing for 
0 and 90 degree is 1mm, while for 45 and 135 degree it amounts to 0 71mm Rays other 
than those from A-scan images are assumed on the basis of symmetry considerations and 
in-between rays are fixed by linear interpolation 
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5.2 EXPERIMENTAL RESULTS 


5.2.1 RESULTS FOR EXCITATION ALONG O" 

In this view foui A-scan images are obtained from experiment performed on LBU set-up 
A-scan plots at the receivers for two representative cases are shown m fig 5 2 Case 1 
and case 2 as shown in fig 5 2, are the A-scan plots at the respective receivers 
corresponding to the pulsed excitations at the extreme left pixel m the 5* and 14*'’ row 

For case 1 the ray (ideal travel path being assumed as linear and coinciding with the 
direction of excitation) has no intercept on the insert region as shown in fig 5 1(c), so the 
TOF obtained from the plot can taken as the time to travel 30mm of glass/epoxy The A- 
scan plot (case 1) shows one sharp signal which initiates at a time of 7 0 ps indicating the 
wave’s arrival 

Thus, the TOF, in this case is taken as 7 0 ps Case 2 corresponds to the ray traveling 
along the 14"’ row of the grid, which has an intercept on the insert zone In contrast to 
case 1, this plot shows signal initiating little later than 7 0 ps, at 7 2 ps, indicating the 
wave’s arrival Thus, the TOF, in this case is taken as 7 2 ps 

Rest of the two A-scan images for excitation along row S*" and 11* are similar to case 1 
and TOF is taken as 7.0 ps Hence, it can be inferred that defect is present between 1 1"* 
and 17* row 

5.2.2 RESULTS FOR EXCITATION ALONG 90® 

In this view four A-scan images are obtained from experiment performed on LBU set-up 
A-scan plots at the receivers for two representative cases are shown in fig 5 3 Case 1 
and case 2 as shown m fig 5 3, show the A-scan plots at the respective receiver ends 
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corresponding to the pulsed excitations at the extreme left pixel in the 5**’ and 14* 
column 

For case 1 the ray (ideal travel path being assumed as linear and coinciding with the 
direction of excitation) has no inteicept on the insert region as shown in fig 5 1(d), so the 
TOF obtained from the plot can taken as the time to travel 30mm of glass/epoxy The A- 
scan plot (case 1) shows one sharp peak which initiates at a time of 7 0 ps indicating the 
wave’s arrival 

Thus, the TOF, in this case is taken as 7 0 ps Case 2 corresponds to the ray traveling 
along the 14* column of the grid, which has an intercept on the insert zone In contrast to 
case 1, this plot shows signal initiating little later than 7 0 ps, at 7 2 ps, indicating the 
wave’s arrival Thus, the TOF, m this case is taken as 7 2 ps 

Rest of the two A-scan images for excitation along column 8* and 1 1* are similar to case 
1 and TOF is taken as 7 0 ps Hence, it can be inferred that insert is present between 11* 
and 17* column Same TOF for both 0° and 90° views signify good experimental results 
as the composite is cross ply and hence nearly same results were expected 

5.2.3 RESULTS FOR EXCITATION ALONG 45® 

In this view nine A-scan images are obtained from experiment performed on LBU set-up 
A-scan plots at the receivers for the nine representative cases are shown m fig 5 4 It can 
be inferred from 0 and 90 degree excitations that the insert is within a 6mmX6mm central 
region of the domain under consideration, so it is expected that only in the last two cases, 
VIZ 8* and 9* case, the ray has an intercept with the insert 

Case 1 to case 7 denotes the signal at receiver end corresponding to pulsed excitation 
from 26* to 8* column on X-axis at a regular spacing of 3 columns between each 
excitation The TOF for each ray is the time of initiation of the signal, signifying the 
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arrival of ray Case 8th and 9* correspond to the rays originating at 2"^' and 5* column 
on the X-axis 

5.2.4 RESULTS FOR EXCITATION ALONG 135® 

TOP data foi 135® aie assumed on the basis of experimental results along 45® The data 
set IS taken same as for 45® 

5.3 RECONSTRUCTION RESULTS 

Reconstruction of slowness field from the TOP data set can be performed by using 
different variations of ART and MART Simple ART is used for reconstruction purpose 
The total number of rays used for reconstruction, as discussed above amounts to 180, 30 
each for 0 and 90 degrees, while 60 each for 45 and 135 degrees In the following 
sections, the effect of change of relaxation parameter, initial guess and algorithm on the 
reconstruction is noted 

The drawback of using isotropic algorithm to study anisotropic material like composite, 
IS to compromise on the actual field values of the material crossection and to some extent 
have to worry about the presence of unexpected defects, but the ease of working with 
simple assumptions such as straight ray approximation, and the successful determination 
of defect location and geometry, along with the anisotropy present makes the study worth 
pursuing Although, some more analysis is required with the assumptions, but simple 
assumption with detailed analysis can be woven into a successful, and more important a 
general tool 
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5.3.1 RECONSTRUCTION USING ART WITH DIFFERENT 
RELAXATION PARAMETERS 


Reconstruction results for diffeient relaxation parameters are shown in figs 5 5(a)-5 5(h) 
For these reconstructions initial guess (/■') is a square matrix with 0 21 as the value of 
each element and stopping criteria is given as, 

abs[(/‘*'-/‘VO]x loose (5 1 ) 

where /**' and /* are the field values after k+1 and kth iteration respectively, 
e IS a suitable stopping criteria, here if s value is 0 1 

Fig 5 5(a) - 5 5(h) are the reconstructions for relaxation parameters ranging from 1 99 to 
0 0001 As can be seen from the reconstructions, the results significantly improve on 
decreasing the value of relaxation parameter and confidence for the presence of central 
circular defect increases on going for values lower than 0 5 Fig 5 5(f) is the 
reconstruction for X=0 01, field values in the central region are m striking contrast with 
rest depicting presence of an extended central defect Figs 5 5(g) and 5 5(h) are for 
X=0 001 and A,=0 0001, a near circular defect of size from 4 to 5mm can be inferred from 
them 

Longitudinal wave velocity m glass epoxy is around 4200m/s and longitudinal wave 
velocity of steel is around 5000m/s (For subsequent analysis these values are taken as 
reference) Figure 5 6(c), depicts some regions where longitudinal wave velocity m glass 
epoxy is exactly matching with nearly 1 8% error in reconstruction for steel 

ART have many variations such as Gordon ART, Gilbert ART etc, difference between 
them lies in the correction value and/or implementation stage Until now, simple ART 
was used for reconstruction, now implementation of Gordon ART and Gilbert ART will 
help in comparing the results Fig 5 6(a)-5 6(d) shows the reconstruction using Gordon 
ART Figure 5 6(d), shows some regions where error m Longitudinal wave velocity m 
glass epoxy is below 10% and nearly 9% error m reconstraction for steel Fig 5 7(a)- 
5 7(d) shows the reconstruction using Gilbert ART On considering Fig 5 7(d) it can be 
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seen that the erroi in longitudinal wave velocity in glass epoxy region is higher than for 
Gordon AR f, while there is a slight reduction in error for insert region reconstruction 
values The corresponding reconstructions using Simple ART figs 5 5(a)-5 5(h), shows 
results nearer to Gilbert ART 

It can be seen that the relative values of longitudinal velocities in the defect and base 
material are on the wrong side, but because of the nature of defect and the data 
acquisition and reconstruction method, defect is considered, as including, also the region 
considerably influenced region with the implantation of insert 

So, It can be said that the TOP data is an information about the influence of affected 
region on the incident ray In that sense some meaningful results can be generated 
without loss of generality Here broken fibres, resin rich zone and bypassing of the ray 
can be attributed to greater slowness values of the defect 

5.3.2 RECONSTRUCTION USING MART WITH DIFFERENT 
RELAXATION PARAMETERS 

Results are here augmented using MART algorithms and observe differences, if any, in 
reconstruction Figures 5 8(a)-(h) are the reconstructions using MARTI for X’s, varying 
from 1 99 to 0 0001 As can be seen from the figs the reconstruction (On considering 
only the defect geometry, condoning the slowness values) didn’t deteriorate for higher 
values of X as obseived using ART and good defect reconstruction is obtained for X as 
high as 1.99, obviously (see below), at such high value of A., it is highly sensitive to imtial 
guess, so our choice of initial guess becomes critical to get good reconstruction for such 
high values, here our initial guess is a square matrix with 0 23 as the value of each 
element Figures 5 9(a)-(h) are the corresponding reconstructions using MART 2 with 
initial guess as 0 21 and nearly similar results are observed 
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Figures 5 10(a)-(f) are the lesults using MART 3 The algorithm has failed to converge 
for X>1.0, while test of the reconstructions are similar to other MART algorithms 

5.4 RECONSTRUCTION WITH DIFFERENT INITIAL FIELD 
VALUES 

Reconstruction results for different initial field values are shown in figs 5 1 l(a)-5 1 1(d) 
For each of these reconstructions Simple ART has been used with 1 fixed at 0 01 and e at 
0 1 ART IS an iterative technique and numbei of iterations to converge depends on initial 
field values, but final field should be neaily independent of initial guess, on using a fixed 
stopping criteria This is clearly evident from the reconstruction results On varying ini tial 
guess field from 0 to 1 0, it is found that the reconstruction is nearly same for all of them, 
except for fig 5 11(c) where / '=(9 The number of iterations to converge for initial guess 
as zero are considerably lower than with otlier cases, and this can be the reason for the 
exception Table 5.1 lists the number of iterations to converge for different initial field 
values 


S.No 

Number of 

itciations (n) 

Initial guess (fi) 

1 

424 

0.2 

2 

232 

1 

3 

32 

0 

4 

176 



Table 5.1 . Depicts Change in the number of iterations to converge 
for different initial field values using simple ART 


The combined effect of relaxation parameter and imtial guess has been studied 
extensively, and the reconstruction is found to be nearly independent of initial guess for 
different values of relaxation parameter This is clearly evident from figs 5 12(a)-(b) that 
as the value of X changes the independence is maintained as discussed above 
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The results using MART aie m striking contrast to that from ART as the relationship 
between icconstruclion and initial field value was affected by the value of relaxation 
paiametcr used figures 5 13(a)-(d) are the reconstructions using MARTI It can be seen, 
that on changing X liom 0 1 to 0 5, the reconstruction is no longer independent of initial 
guess Similai lesulls using MART 2 are shown in figs 5 14(a)-(d) 

5.5 NUMERICAL SIMULATION 

To check the performance of algoiithms used above, the dataset is generated using a 
simulated specimen similar to the expeiimental object, i e a 3 0X3 0mm specimen with a 
circular central insert of diameter 6mm The reconstruction is carried out by considering 
weight matrix obtained from stiaight-lme assumption The slowness of the insert is 
taken as twice of base material The longitudinal velocity in the defect free region is 
taken as 30mm/}.LS The no. of views as well as the no of rays per view remain the same as 
used for experimental results. To be able to compare better with the experimental results 
the error levels are calculated for all the algorithms with different values of relaxation 
parameter The results aie presented under the following subheadings The error ELI and 
rras error EL2 for subsequent error analysis are as follows, 
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5.5 1 Numerical Simulation using ART Algorithms 


To be able to compare better with the experimental results the error levels are calculated 
for all algorithms used for different values of relaxation parameters 

1. Simple ART 


SNo 

ELI 

EL2 

1 

Iterations 

1 

11 32 

16 53 

1 99 

401 

2 

11 33 

14 8515 

1 5 

71 

3 

11 37 

16 9 

1 0 

313 

4 

10 755 

1621 

05 

114 

5 

10 56 

15 99 

0 1 

477 

6 

14 56 

27 08 

0 01 

111 

7 

14 56 

25 8 

0 001 

1198 

8 

661 

1274 3 

0 0001 

224 


Table-5 2 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using Simple ART 


Different error levels for a range of relaxation parameters are presented for Simple ART 
The algorithm doesn’t seem to follow any particular dependence of error or number of 
Iterations on X . Minimum error levels are observed at =0 1 
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2. Gordon ART 


SNo 

ELI 

EL2 

A 

Iterations 

1 

8 19 

10 74 

1 99 

54 

2 

1061 

15 76 

1 5 

61 

3 

10 56 

15 72 

10 

74 

4 

10 61 

15 786 

05 

570 

5 

9 71 

14 85 

0 1 

186 

6 

12 99 

22 04 

001 

129 

7 

20 82 

47 06 

0 001 

465 

8 

91 61 

361 82 

0 0001 

566 


Table-5 3 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using Gordon ART 

Different error levels for a range of relaxation parameters are presented in table 5 3 for 
Gordon ART It has shown better performance than Simple ART for all values of X, with 
significant reduction in error levels The algonthm is most efficient in terms of number of 
iterations and error values for X=0 1 

3. Gilbert ART 


SNo 

ELI 

EL2 


Iterations 

1 

85 62 

87 23 

05 

4 

2 

12 23 

19 9 

0 1 

15 

3 

12 16 

19 86 

0 01 

158 

4 

29 82 

80 92 

0 001 

258 


TabIe-5.4 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using Gilbert ART 
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Gilbert ART converged the problem only for relaxation parameters lesser than 0 5 and 
except for X,=0 1 and A.=0 01, the error values are too high 

5.5 2 Numerical Simulation using MART Algorithms 

1. MARTI 


SNo 

ELI 

EL2 

X 

Iterations 

1 

0 9083 

8 9665 

1 99 

5 

2 

21 59 

27 57 

1 5 

3 

3 

39 97 

41 54 

05 

3 

4 

9 56 

1421 

0 1 

15 

5 

10 41 

15 876 

0 01 

118 

6 

26 159 

77 808 

0 001 

169 

7 

68 68 

353 33 

0 0001 

381 


Table-5 5 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using MART 1 

MART 1 IS different from the above ART algorithms in the sense that it has shown 
least error levels at such a high value of relaxation parameter with few iterations 
Also for A, =0 1 the performance is better than ART algorithms 
2. MART 2 


SNo 

ELI 

EL2 


Iterations 


28 42 

31 33 

1 99 

4 


35 29 

37 94 

1 5 

9 

^Hi 

87 93 

88 62 

05 

17 

am 

8 42 

13 42 

0 1 

36 

5 

11 28 

19 93 

0 01 

136 

6 

50 29 

172 83 

0 001 

138 

1 

665 56 

3221 3 

0 0001 

77 


Table-5.6 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using MART 2 


53 




Different error levels for a range of relaxation parameters are presented in Table 5 6 foi 
MART 2 The algorithm is most efficient m terms of number of iterations and error 
values for A,=0 1 and also has the least error level for X=0 1 among the additive and 
multiplicative category 

3. MARTS 


SNo 

ELI 

EL2 

2 

Iterations 

1 

9 25 

13 66 

0 1 

18 

2 

13 17 

16 06 

001 

40 

3 

1795 

21 81 

0 001 

140 

4 

102 63 

404 12 

0 0001 

79 


Table-5 7 Different error levels, number of iterations, and relaxation parameters are given 
for the reconstruction of simulated specimen using MART 3 

MART 3 failed to converge for relaxation values more than 0 1 and have shown best rms 
EL2 and ELI for 1 

Errors due to reconstructions for these six algorithms are presented The best value of 
relaxation parameter that gives minimum error and fastest convergence has been 
determined through numerical experimentation and has been reported for the ART and 
MART algorithms 

High value of ELI and EL2 for these algorithms can be attributed to limited number of 
views, m this case 4 views Increasing the number of views generally decreases rms error 
and error as low as 2 0 is reported for 8 views It is remained as such to be able to 
compare better with the experimental results Gilbert ART have shown largest RMS 
error In ART family, Gordon ART has shown best results for this problem 
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MART algorithms have shown less error (ELI as well as EL2) than ART family for 0 01 
< A- < 0 1 Particularly, MART 2 has shown best peiformance at A=0 1 m the additive and 
multiplicative family of algorithms 

5.6 SIMULATED RESULTS 


25 

20 

1S 

10 

5 

5 10 15 20 25 30 

Simulated Specimen 

Here, simulation is performed to check the performance of the algorithms used Figure 
5 15 to 5 20 are the results obtained after reconstruction The reconstructions using 
MARTI and MART2 are better than other algorithms as expected from their RMS 
values, so these two algorithms are shortlisted for further analysis 



5.7 REVIEW OF EXPERIMENTAL RESULTS BASED ON 
SIMULATED STUDY 

The superimposition of experimental results with the simulated study clears the factor of 
algorithm dependency on confusion of image interpretation, so that physical reasons can 
be ascertained to deviations from the examined body The experimental results are again 
discussed to fix some of the suspicions MARTI and MART2 are used for their better 
performance as discussed above for = 0 1 
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The lines at 45 and 135 degrees in the experimental result can now be ascertained to the 
directional dependence or anisotropy m the base material, as the result of application of 
anisotropic algorithm on an anisotropic material Similar results for 0 and 90° affirms the 
cross ply arrangement of fabric 

As can be seen from the results (Fig 5 21 to fig 5 24) that it can generate defect location 
and geometry The size of the defect is somewhat large, though not far away from the 
actual values The reason can be attributed to the bending of ray near the edge of the 
flaw, which IS not accounted for m the straight line assumption 

5.8 POSSIBLE SOURCES OF ERROR 

The data acquired from laser based Ultrasonics (LBU) setup are for limited views and 
few number of rays per view, so the following sources of errors are inevitable 

1 Experimental error 

• Location 

Fine manual adjustments are required to fix the position of sources and 
receivers to the marked points and the inaccuracies in positioning may 
result m error in subsequent reconstruction 

• Time Of Flight (TOF) 

There may be some error in taking the Time Of Flight (TOF) value from 
the A-scan images 

2 Incomplete data 

Linear interpolation has been used to fix the m between projections and this 
may be one of the sources of error 


56 



5.9 Closure 


In this chapter, the Time Of Flight (TOF) information from the A- Scan images has been 
obtained from Laser-Based Ultrasonic (LBU) set-up and subsequently reconstructed the 
slowness values for the domain under consideration using different ART and MART 
algorithms for different values of relaxation parameter and initial guesses and compared 
the results with numerical simulation 
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Case 1 - Input at 5* row 



Case 2 - Input at 14* row 


Figure 5 2 Signal at the receiver end due to excitation along 0® (X-direction) 
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Case - 3 Time of Flight (4 0 ^s) 






Case - 4 Time of Flight (7 4 p.s) 



Case - 5 Time of Flight (9 0 (is) 
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Case - 9 Time of Flight (17 5 |j.s) 

Figure 5 4 Signal at the receiver end due to excitation along 45° wrt X-axis 






Figure 5 5(e)- A,=0.1 


Figure 5 5(f); >,=0.01 



Figure 5.5(g); >,=0.001 Figure 5.5(h); >,=0.0001 


Figure 5.5(a)-5.5(h): Reconstruction of slowness values for different relaxation parameters 
using Simple ART (Material: Glass/epoxy with central steel insert) 



Figure 5.7(a) A,=0.1 


Figure 5.7(b) X=0 01 



Figure 5 7(c) A,=0 001 


Figure 5.7(d) X=0.0001 


Figure 5.7(a)-5.7(d). Reconstruction of slowness values for different relaxation parameters 
using Gilbert ART (Material. Glass/epoxy with central steel insert) 










Figure 5.10(e): ^=0 001 


Figure 5 10(f) A,=0 0001 


Figure 5.10(a)-5.10(f): Reconstruction of slowness values for different relaxation 
parameters using MART 3 with =0.21 (Material: Glass/epoxy with central steel insert) 











Figure 5.12(a). /=0 1 
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Figure 5.12(b) /'=0.21 


Figure 5.12(a)-5.12(b): Reconstruction of slowness values for different initial guess 
(A, =0.5, IS same for all the cases) Material is Glass/epoxy with central steel insert. 


Figure 5 13(a) /=0 21, A,=0 1 


Figure 5 13(b) / =0.1, A,=0 1 



Figure 5.13(c)- /=0 1, A,=0.5 Figure 5.13(d) :/=0.21, X,=0.5 


Figure 5.13(a)-(d) Depicting the combined effect of initial guess and relaxation parameter 

on reconstruction, using MART 1. 
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Figure 5.14(b):/=0.21, X=0 5 
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Figure 5.14(c) :/'=0.1, A.=0.1 
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Figure 5.14(d)* /=0 21, A.=0.1 



Figure 5.14(a)-(d) Depicting the combined effect of initial guess and relaxation parameter 

on reconstruction, using MART 2. 


Fig 5 15 Reconstiuction using Simple ART 
(30 rays, 4 views, X = 0 1) 


Fig 5 16 Reconstruction using Gordon ART 
(30 rays, 4 views, Jl = 0 \ ) 


Fig 5 17 Reconstruction using Gibert ART 
( 30 rays, 4 views, >1 = 01) 


Fig 5 18 Reconstruction using MARTI 
(30 rays, 4 views, >1 = 01) 
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Fig 5 19 Reconstruction using MART2 
(30 rays, 4 views, >1 = 01) 


Fig 5 20 Reconstruction using MART3 
(30 rays, 4 views, >1 = 01) 







Chapter 6 

CONCLUSIONS AND SCOPE FOR FUTURE WORK 


6.1 CONCLUSIONS 

In the present study of tomographic reconstruction of defects in composite materials, 
most of the ART and MART algorithms are used for reconstruction and thereafter the 
results are compared with simulated reconstructions Laser based ultrasonic setup is used 
to find the time of flight data Based on the results obtained, the following conclusions 
are drawn 

(1) Carefully obtained Time Of Flight (TOF) data from Laser based ultrasonic 
(LBU) setup IS sufficient and satisfactory enough to be used for further Non- 
Destructive analysis 

(2) The Algorithms (ART and MART) used here, neglecting ray bending, is capable 
to generate defect location and geometry The size of the defect may come 
somewhat larger than the actual value 

(3) The reconstruction is found to be nearly independent of the initial guess value, 
which indicates the considerable insensitivity of the Algorithms to initial 
approximation, except for MART, when high values of relaxation parameter 
were used 

(4) MART algorithms, especially MART 2 shows better performance than additive 
algorithms 

(5) All the additive and multiplicative algorithms we discussed had shown less rms 
error for X=0 1, hence making it the preferred relaxation parameter for these 
algorithms 

(6) The superimposition of simulation study on experimental results can enhance our 
confidence over the anisotropy in object 
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(7) Finally, the straight ray approximation and linear-interpolation, being 
oversimplifications, have not interfered with the accuracy of reconstruction to an 
extent, which may demand necessary incorporation of complexities 

6.2 SCOPE FOR FUTURE WORK 

(1) The time of flight data obtained from shear waves can be used for finding the 
defects in composite materials 

(2) Experimental determination of shear waves is possible by making some changes 
m the interferometry 

(3) Non-linear interpolation can be used to complete the data set 

(4) Algorithms such as Maximum Entropy and Mimmum Energy can be tried for 
reconstruction 
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Avvendix-A 

SIMPLE ART PROGRAM USING MATLAB 


% art z IS convergence factor 
%s IS 30*30 slowness matrix 
% err is stopping criterion 

% diff is difference between projection and calculated data 

% y is number of iterations 

diff=0,z=0 001,err=l 0,y=0, 

s=ones(30,30),s=0 21 *s, 

s0=s, 

% for 0 deg 
while(err>=0 1), 
s0=s, 

fid=fbpen('dataaO txf ,'rt'), 

a=fscanf(fid,'%f ),i=l , 

while(i<=30), 

sum=0, 

forj=l 1 30, 

sum=s(i,j )* 1 +sum, 

end 

diff=a(i)-sum, 
forj=l 1 30, 

s(ij)=s(i,j) + z*l*(diff/30 0), 

end 

1 = 1 + 1 , 

end 

fclose(fid), 

% for 90 deg 

fidl=fopen('data90 txtVr'), 

a=fscanf(fidl,'%f), 

diff=0,j=l, 

while(j<=30), 

sum=0, 

for 1=1 1 30, 

sum= 1 * s(i j)+sum, 

end 

diff=a(])-sum, 
for 1=1 1 30, 

s(ij)=s(ij) +z*l*(diff/30 0), 

end 

J=J+1, 

end 

fclose(fidl). 
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% for 45 and 135 deg 

fidl=fopen('data45 txtVr'), 

a=fscanf(fidl,'%f), 

diff=0,diffl=0, 

for k=l 1 30, 

i=k, sum=0,w=0,suml=0, 

forj=l 1 k, 

sum=sum+0 707*s(i,j),suml=suml+0 707*s(31-i,j), 
w=w+0 707, 
ifi>l 
1 = 1 - 1 , 

sum=sum-t-0 707*s(i,j),sunil=suml+0 707*s(31-i,j), 
w=w+0 707, 

end 

end 

diff=a(k)-sum, 

diffl=a(k)-suml, 

i=k, 

forj=l 1 k, 

s(ij)=s(ij) +z*0 707*diff/w, 
s(31-i,j)=s(3 l-ij)+z*0 707*diffl/w, 
ifi>l 
1 = 1 - 1 , 

s(ij)=s(i,j)+z*0 707*diff/w, 
s(3 l-i,j)=s(3 l-ij)+z*0 707*diffl/w, 

end 

end 

end 

fork=l 1-30, 

j=k, smn=0,w=0,suml=0, 

for 1=30 -l*k, 

siim=sum+0 707*s(ij), 

sum 1 =sum 1 +0 707 * s(3 1 -i j ), 

w=w+0 707, 

ifj<30 

J=J+1> 

sum=sum+0 707*s(ij), 
suml=sumH-0 707*s(31-i,j), 
w=w+0 707, 

end 

end 

diff=a(k+3 0)-sum, 
diffl=a(k+30)-suml , 

J=k, 
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for 1=30 -Ik, 

s(i,j)=s(i,j)+z*0 707*diff/w, 
s(31-i,j)=s(31-ij)-i-z*0 707*diffl/w, 

ifj<30 

s(i,j)=s(ij)+z*0 707*diff/w, 
s(31-i,j)=s(31-i,j)+z*0 707*diffl/w, 

end 

end 

end 

eir=abs((sO-s)/(s)* 1 00), 
y=y+l, 

end 

% change of units to m/sec 
s=s*0 001, 

% tomographic image generation 

r=l 1 30, 

z=s, 

[x,y]=meshgrid(r,r), 

colormap(hsv), 

pcolor(x,y,z), 

rotate3d on, 

view(3), 

grid, 

shading mterp, 
colorbar, 
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